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Abstract 

r~| ■ Certain dark matter interactions with nuclei are mediated possibly by a scalar or pseudoscalar 

Higgs boson. The estimation of the corresponding cross sections requires a correct evaluation of 
the couplings between the scalar or pseudoscalar Higgs boson and the nucleons. Progress has been 
^ . made in two aspects relevant to this study in the past few years. First, recent lattice calculations 

show that the strange-quark sigma term 0"^ and the strange-quark content in the nucleon are 
much smaller than what are expected previously. Second, lattice and model analyses imply sizable 
SU(3) breaking effects in the determination on the axial-vector coupling constant g\ that in turn 
affect the extraction of the isosinglet coupling g\ and the strange quark spin component As from 
polarized deep inelastic scattering experiments. Based on these new developments, we re-evaluate 
the relevant nucleon matrix elements and compute the scalar and pseudoscalar couplings of the 
proton and neutron. We also find that the strange quark contribution in both types of couplings 
is smaller than previously thought. 
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I. INTRODUCTION 



Knowledge of scalar and pseudoscalar Higgs boson interactions with the nucleons at low energies 
is an important ingredient for computing certain dark matter-nuclei interaction cross sections. For 
example, the scalar interactions between fermionic dark matter and nucleons mediated by the scalar 
Higgs boson can be studied by measuring spin-independent (SI) cross sections, while interactions of 
the dark matter with the nucleons mediated by light pseudoscalars can be probed by studying the 
momentum-dependent, spin-dependent (SD) cross sections. ^ At the quark level, the interactions 
of a neutral scalar cj) and a pseudoscalar a with quarks have the general expressions: 

q 

C^gg = (^/2Gi.)i/2^^gmggn5#, (1) 
q 

where and are the couplings of the scalar (j) and the pseudoscalar a with quarks, respectively. 
While in the standard model = 1 for all g's, they may have values different from unity beyond 
the standard model. To evaluate the effective scalar and pseudoscalar couplings with the nucleons 
given by 

94NN = (V2Gir)i/2^(Ar|c,mgg(?|iV), and 
q 

9aNN = {V2GFf'^Y.^N\iqmqqi^^q\N), (2) 
q 

it amounts to computing the nucleon matrix elements of quark scalar and pseudoscalar densities, 
namely, {N\mqqq\N) and {N\'mqqi'y^q\N) . 

The coupling q^nn was first studied by Shifman, Vainstein and Zakharov (SVZ) [l| who assumed 
a negligible strange quark contribution to the nucleon mass. Consequently, their scalar-nucleon 
coupling was dominated by heavy quarks. Based on the pion-nucleon sigma term available in late 
80's which implied a sizable strange quark content in the nucleon, it was shown in 0, S] that the 
effective couphng g^^ was dominated by the . quark rather than by heavy quarks. As a result, 
the scalar-nucleon coupling was enhanced by a factor of about 2.5 [3|. 

Although the coupling of the axion, which is an example of the pseudoscalar boson, with the 
nucleon had been studied before (see Q] for a review), a thorough discussion of generic pseudoscalar- 
nucleon couplings was only given Me7 in a based on .he relation 

{N\ui-i^u + di-i^d + si-i^s\N) = , (3) 

derived from the large- A^c and chiral limits. Contrary to the case of scalar-nucleon couplings, g^NN 
does receive significant contributions from light quarks even in the chiral limit. It is conventional 
to express gaNN in terms of the axial-vector couplings g\ {a = 0, 3, 8) or the quark spin components 
Ag with q = u,d, s. 



Since the cross sections for the dark matter interacting with the nucleons induced by pseudoscalars vanish 
in the non-relativistic limit, it means that the signature of dark matter cannot be seen by virtue of SD 
direct detection. 
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Recently, there is some progress in topics relevant to the interactions between the scalar or 
pseudoscalar Higgs boson and the nucleons. First, there have been intensive lattice calculations of 
the pion-nucleon sigma term a-j^Ni the quark sigma term Uq = mq{p\qq\p) , and the strange quark 
content in the nucleons characterized by the parameter y. Especially, the lattice calculation of ag 
is not available until recent years. The lattice results indicate that the strange quark fraction in 
the nucleon and Cs are getting smaller than what we had one or two decades ago. For example, 
the new lattice average of = (43 it 8) MeV 0] is much smaller than the value of ~ 390 MeV in 
the late 80's and ~ 130 MeV in the early 90's. Second, recent lattice and model calculations 0] 
hint at a size of g\ about 20% smaller than the canonical value 0.585 determined from the hyperon 
/3-decays supplemented by flavor SU(3) symmetry. Since g\ is extracted from the measurement 
of the first moment of the proton polarized structure function g\, a decrease in will lead to an 
increase in g\. This in turn implies a As reduced in magnitude by a factor of 2 to 3. 

Motivated by the above-mentioned developments, in this work we would like to pay a revisit to 
the couplings of scalar and pseudoscalar bosons with the nucleons to incorporate the recent progress. 
Moreover, we wish to investigate if it is possible to evaluate the matrix elements (A^|GG|iV) and 
{N\qi'y^q\N) without invoking the large- A'^c chiral relation Eq. ([3]) as this relation is presumably 
subject to l/Nc and chiral corrections. 

This paper is organized as follows. Section |Tl] deals with the couplings between the nucleons 
and the scalar Higgs boson, whereas Section |TTI] analyzes the case with a pseudoscalar Higgs boson. 
Our findings are summarized in Section IIVI Appendix [SI provides a comparison of our results 
against two sets of parameters commonly used in dark matter physics analyses. 



II. SCALAR HIGGS COUPLINGS TO THE NUCLEONS 



The baryon mass can be expressed in terms of the matrix element of the trace of the energy- 
momentum tensor 

9^^ = niuuu + niddd + nigSS + ^ mhqhqh + -T—^G , (4) 

h ^"'^ 

with 

5q,2 2 2 

/3 = h = ll- -n^ - -nh , (5) 

where n^^/i denote the numbers of light and heavy quarks, respectively, and the trace anomaly is 
governed by the GG = G'^^G"'^'^ term. As pointed out by Shifman, Vainstein and Zakharov [l|, 
heavy quarks contribute to the baryon mass via the triangle diagram with external gluons. Under 
the heavy quark expansion [l| 

„,j,,,^-^^GG + o(i|) , (6) 

where is a typical hadron mass scale. Thus we have 

Q^ii^ = niuUU + niddd + nisss — - — GG. (7) 
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Under the light quark mass expansion, the baryon mass has the general expression 

<? g q 

3/2 

where the coefficients Bg, Cg and Dg are baryon state-dependent. The rug terms arise from chiral 

^ 3/2 

loop contributions which are ultraviolet finite and non-analytic in the quark masses as nip oc nig 



with mp being the pseudoscalar meson mass in the loop js-ll]. All terms of order have been 



calculated in |12l . Il3l |. We first consider the quark mass expansion of the baryon masses to the 
leading order. Using the proton matrix elements defined by 

B^ = {p\uu\p), Bd = {p\dd\p), Bs = {p\ss\p) (9) 

and applying the fiavor SU(3) symmetry, the masses of the octet baryons to the first order can be 
expressed as 

nip = nio + Bumu + Bdnid + Bgnig, 

tnn = niQ + Bdniu + Bunid + BsiUs, 

mj^+ = niQ + Bumu + Bsnid + Bdms, 

ni^o = mo + {Bu + Bs)ni + Bdnis, (10) 

ms- = nio + Bsniu + Bunid + BdiUs, 

m^o = niQ + Bdniu + Bsnid + Bums, 

m^- = mo + Bsmu + Bdnid + BuiUs, 

1 1 

mA = mo + -{Bu + 4:Bd + S,)?n + -(2B„ - Bd + 2Bs)ms, 

with m = {niu + nid) /2. From the above mass formulae one can derive the Gell-Mann Okubo mass 
relation in the baryon sector 

1 3 

rriE - TTiAT = -(ms - mAf) -(mg - ttia) , (11) 
the Coleman- Glashow mass formula 

nip — nin + nij^- — mj^+ + ni-=o — ni-=- = , (12) 
and the constraint such as 

m=o — ms — 

with TTiAT = {nip + m„)/2, m-£ = (mv+ + ni^-)/2 and ?7t,= = (m=o -|- m=-)/2. 
Using the current quark masses 15 1 

?n„ = 2.5toi MeV , = 5.0t[J;^ MeV , = WOtlo MeV , (14) 

normalized at the scale n = 2 GeV, we obtain ^ 

Bu-B. = ^(^^-^^) ^3.94, 
2ms -mu-md 

Bd-Bs= ^^^--^-^ =2.64 (15) 
2ms - mu-md 



^ From the mass formulae given in ()10p . we also have the relations B„ — Ba = {mY,+ —rn^- )/(mu—rnd) = 3.23 
and Bd—Bs = (m^o — m5;-)/(m„ — md) =2.74. It is clear that the value of _B„ — i3s obtained in this manner 
differs sizably from that in Eq. (fT5|) . Indeed, in this case one needs to take into account electromagnetic 
corrections to the isospin splitting. 
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at ^ = 2 GeV. To determine the parameters B.^, and Bg we need additional information. To 
proceed, we define two quantities: 

which characterizes the strange quark content in the proton and 

do = m{N\uu + dd- 2ss\N) , (17) 
which is related to the pion-nucleon sigma term ^ 

a.„N = rh{N\uu + dd\N) (18) 

by 

(y-KN = ■ (19) 
1 - y 

Hence, ctq is the pion-nucleon sigma term in the absence of the strange quark content in the nucleon. 
It follows from Eq. (llOp that at the leading order ^ 

3m m ^ \ /r.^\ 

(Tq = — [m^ — m\), or (Jq = — (m= + ms — 2mN), (20) 

nis — m rus — m 

Numerically, do = 23.7 MeV and 24.7 MeV, respectively. The experimental and/or theoretical 
information of cTttTv will enable us to determine the strange quark fraction y which puts an additional 
constraint on the parameters B^, B^ and Bs- 

In the early study based on the sigma term Gj^n = 55 ~ 60 MeV, one has a large strange quark 
content y 0.55 and 

an ^ 20 MeV, ad « 33 MeV, a, w 394 MeV, (21) 

where cjg is the scale-independent quark sigma term defined by 

aq = mq{p\qq\p) . (22) 

The strange quark content y and the sigma term ag are both unexpectedly large. 

In recent years, there exist many lattice results on the strange-quark sigma term ag either 
based on a direct calculation 16h20|] or relied on the Feynman-Hellman theorem through a study 



of how the nucleon mass varies with 
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23( 1 . A global fit to the current dynamical lattice 



data yields Ug = 43 it 8 MeV [5|. Meanwhile, the strange quark content also becomes smaller, 
y = 0.05 ~ 0.20. Then the question becomes "Is it possible to reduce the naive estimate of ag by 
one order of magnitude?" This possibility has been realized by several authors [s-lS] by considering 



^ The pion-nucleon sigma term is the scalar form factor cr(t) at t = defined by {N {p')\rh{uu + dd)\N {p)) 
u{p')a{t)u{p). 
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The sigma term ctq can also be expressed as 



nid + m„ 1 
fJo = Toe- - TOs+ - ■7;\^ri - "ipj 



However, it is necessary to take into account the electromagnetic corrections to the isospin splittings 14 1. 
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the corrections to the baryon masses from higher order terms in the hght quark mass expansion. 
Since ah corrections to the baryon masses up to the quadratic order in the hght quark masses have 
been completely computed by Borasoy and Meii3ner based on he avy baryon chiral perturbation 
theory [l^.fist. we shall follow their results. From Eq. (80) of Ref. 13], the quark mass expansions 
of the baryon masses in the fashion of Eq. ([8|) read 



niN = mo(l + 0.34 - 0.35 + 0.24), 
ms = mo(l + 0.81 - 0.70 + 0.44), 



niA 



:mo(l + 0.69 - 0.77 + 0.54), 
?no(l + 1.10- 1.16 + 0.78). 

3/2 



(23) 



Evidently, there is a large cancellation between the BqUiq and non-analytic Cqniq terms. The 
quadratic quark mass contribution is fairly small for the nucleons but sizable for A, S and H. 
Consequently, Eq. ()20p is modified to 

3m 



m 



m 



or 



m 



m 



mo(-1.16 + 0.78) - lUA + mo(-0.77 + 0.54) 
ms + ms - 2mN + mo(1.16 - 0.78 + 0.70 - 0.44 - 0.70 + 0.48) 



(24) 



Using the value of niQ = 767 MeV fitted to the baryon mass spectrum |13[ |. we arrive at ctq = 37 



MeV using either formula in Eq. ()20p . Indeed, the result of ctq = 36 + 7 MeV is quoted in Ref. 13l].^ 
Note that this value is larger than those quoted below Eq. (j20p by about 50%. 

Based on the same approach, Borasoy and Meifiner have also computed the pion-nucleon sigma 
term using the relation aj^N = fa{dmN / dm) and obtained a-j^N = 45 MeV. Substituting this to 
Eq. (fT9]l leads to the strange fraction y = 0.20. Since + o"d w UttTv and mo = 767 MeV, it follows 
that as ~ 126 MeV. Similar results with g^^n = 45 + 7 MeV, ctq = 35 + 5 MeV, y « 0.22 and 
~ 130 MeV were also obtained by Gasser, Leutwyler and Sainio 0]. The pion-nucleon sigma 
term is related to the measured on-shell pion-nucleon scattering amplitude S,riv at the Cheng- 
Dashen point by 

SttAT = CJ,rAr(2m^) + A/j = Cr7rAr(0) + A^ + A/j . 



(25) 



It was found that A^ = 15 MeV [9|] and A/? 2 MeV 

In short, the analyses by Gasser, Leutwyler and Sainio and by Borasoy and Meifiner lead to 
cTttTv ~ 45 MeV, y ~ 0.2 and as ~ 130 MeV. As mentioned before, the recent intensive lattice 
activities in the calculations of a-,^N-, V and in particular as have put previous analyses on test. It is 
obvious that the predicted as is still too large compared to the recent lattice average 43 + 8 MeV. 
If as is fixed to be 50 MeV, then Bs = 0.5 at = 2 GeV. Prom Eqs. dUD and we have 



c^ivN = <7o + 2Bsrh . 



(26) 



Hence, a smaller as means that the difference between a.,^N and a^ should be smaller. With (Tq = 36 
MeV, the parameters (TttTv and y are fixed to be a-j^N = 39.75 MeV and y = 0.094 . Indeed, most 
recent lattice calculations in Refs. [3] and 



23] yield a^N = 38 + 12 MeV and 39 + 41:^7 MeV, 



3 /2 

Because of the strong cancellation between terms of order to, and mq , the resultant cto and (Tjrjv will 
become too large if the corrections to the baryon masses are not included. 
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TABLE I: Summary of parameters for three choices of ag. All the terms are in units of MeV except 
the dimensionless y and g,f,NN- The sigma term gq is fixed to 36 MeV. The extraction method is 
described in the main text. For the neutron, 0"!"^ = 11 MeV and o"^"^ = 32 MeV for (Tg = 50 MeV. 
To get the numerical value of g^N we have assumed three heavy quarks and Q = 1. 





0"m 






{vmGG\p) 


y 


g^NN 


60 


16.2 


21.6 


40.5 


-187 


0.11 


1.2 X 10-3 


50 


16.0 


21.1 


39.8 


-189 


0.09 


1.1 X 10-3 


40 


15.7 


20.6 


39.0 


-192 


0.08 


1.1 X 10-3 



respectively. To conclude, as becomes smaller because of two reasons: a significantly larger ao due 
to higher-order quark mass corrections to the baryon masses and a smaller Gt^n and its difference 
with (To. 

To determine the parameters and B^, we apply the higher-order quark mass corrections 
given in Eq. (|23p to Eq. (|15p and get 

Bu-Bs ^ ms - mjv + mo(1.16 - 0.78 - 0.35 + 0.24) ^ ^ 

Bd-Bs ms -mjv + mo(0.70 -0.44 - 0.35 + 0.24) ' ' ^ ' 

This together with the sigma term (T,rAf = 39.75 MeV yields (T^ = 16.0 MeV and cr^ = 21.1 MeV. It 
should be stressed that these parameters are all correlated. For example, one cannot have a small 
as of order 50 MeV and a (TttAt larger than 40 MeV. Various parameters are summarized in Table U] 
with as varies from 40 MeV to 60 MeV. 

We are now ready to compute the effective scalar-proton coupling strength 



= (\/2Gi7')^/^ {^u(yu + CdCTd + Cs(^s + -^{rrip - au - ad - as) ^ Chj , (28) 

where use of Eqs. ^ and ([7]) has been made. If light quark contributions are neglected, one will 
have 

g^N ^ {V2GFf'^^ niN Y <h , (29) 

h 

- i 

which is the original SVZ result [Ij]. Assuming three heavy quarks and = li we then have 
gcfipp = 1.1 X 10 3 from Eq. (128p which is smaller than the previous estimate 2.2 x 10-^ by a 
factor of 2. Note that since the difference between g(f,nn and gi^pp is negligible, they can be denoted 
collectively as g^NN- We also note in passing that the value of g^NN depends on the number of 
heavy quarks. The increment is 2.6 x 10-^ for each additional heavy quark. 
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III. PSEUDOSCALAR HIGGS COUPLINGS TO THE NUCLEONS 



The analog of the heavy quark expansion ([6]) in the pseudoscalar sector is [l| 

mhqhh^qh ^ -^GG + O (^-^^ , (30) 
with GG = ^e^i/a/jG^^'^G""'^. The pseudoscalar coupling with the nucleon then becomes 



9aNN 



^GG 
Svr 



N) \ , (31) 



where Eg = {N\qi^^q\N) . The quantities Eg are related to the nucleon matrix elements of axial- 
vector currents. The divergences of the nucleon matrix elements of the axial-vector currents 



A,, = U7fi75u + d-f^-y5d - 257^755 , (32) 

read (we have dropped i^i^^tp for convenience) 
Qa^n = ni-uEu- uidEd , 

g\mM = m^E^ + m,E, - Im^Es , (33) 



g^rriN = rriuEu + maEd + mgEs + 3 ( iV 



^GG 
Svr 



N 



and g'\ are the axial-vector form factors of {N\A'^\N) at q^ = 0. Owing to the gluonic anomaly, it 
is clear from Eq. (|33|) that E^, E^ and Eg cannot be determined individually without additional 
information on the matrix element (A^|GG|A^). In the OPE approach, the proton matrix element 
of the axial- vector current is related to the quark spin component: ^ 

{p\Qliil5q\p)s^ = Ag , (34) 

where is the proton spin 4- vector and Aq represents the net helicity of the quark flavor q along 
the direction of the proton spin in the infinite momentum frame 

Aq= [ Aq{x)dx = [ Iq^ix) + q^{x) - q^{x) - q^ix)] dx . (35) 
The axial coupling constants can be expressed as 

g\{Q^) = Au{Q^) - Ad{Q^) , 

g\{Q^) = Au{Q^) + Ad{Q^) - 2As{Q^) , (36) 
g'iiQ^) = Au{Q^) + Ad{Q^) + As{Q^) , 



^ In the late 80's and early 90's there had been a hot debate whether or not there was an anomalous gluon 
contribution to (p|97^759|p) and to the first moment of the polarized structure function g\ through axial 
anomaly. It was realized later that it depended on the factorization scheme in defining the quark spin 



density Aq{x) (for a review, see 2J1). Eq. (p4)) is valid in the gauge- invariant factorization scheme such 
as the MS scheme. 
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with Q2 = _^2_ wj^iig the non-singlet couphngs and are readily determined from low- 
energy neutron and hyperon beta decays, the singlet coupling can be inferred from experiments 
on polarized deep inelastic lepton-nucleon scattering. The measured first moment of the polarized 
proton structure function gi{x) is related to g'\ via 

r?(Q2) ^ £ g{{x, Q')dx = Cns(Q') [^9\ + ^b'a) + ^Cs(Q^)5^(Q2) , (37) 

where Cns and Cs are the perturbative QCD corrections to non-singlet and singlet parts of r^*, 
respectively. They have been calculated to order a^, and the explicit expressions can be found in 
25|. 



Historically, there had been several attempts for solving Eq^s from Eq. (|33p : 

1. Inspired by the early EMC result of a small g^ 26[, it was natural to consider the large- A^c 



chiral limit (i.e., the double limits Nc ^ oo and mq — )• 0) (27|, l28| where 5^4 = and 



Eu + Ed + Es = . (38) 

However, later deep inelastic scattering experiments showed that g^ was of order 0.30 or 
even bigger. This means that the relations ~ and Eq. (|38|) are subject to 1/Nc and 
chiral corrections. 



2. Some authors, e.g., 3l|, IS^], in early 90's had tried to relate {p\GG\p) with the gluon spin 



component AG of the proton; more precisely, 

2mNAGui-f5U = {p\GG\p) . (39) 

A measurement of AG will then allow us to extract Eg^s. However, the above relation is 
repudiated on the ground that there is no twist-2 gauge-invariant local operator definition 
for AG. Indeed, if we write = — UfK^ + UfK^ = + UfK^ and notice that 
^^,K^' = {as/2TT)GG and 

{N{p')\K^\N{p)) = n(p')[Gi(g2)7^75 + G2{q^)q^7Mp), (40) 

we shall see that the forward matrix element of d^K^ does receive a contribution from 
G2{<f')q^ even in the limit — )• as a solution of the U(l) problem demands a ghost pole 
coupled to K^^. Its value depends on the way of taking the forward limit — )• so that it is 
indefinite. Another serious problem is that Eq. (I39p will imply a large isospin violation for 
the gluon spin component AG; that is, (AG)n is very different from (AG)p since (A^|GG|iV) 



violates isospin symmetry [33|, l3j] . 



3. It is known that in the presence of strong GPviolation, the GP-odd operator —(g^ /32-k'^)6GG 
in QCD can be rotated away in terms of quark pseudoscalar densities, i.e., the well-known 



Baluni operator |35| | 



^^Baiuni ^ Qjn^ui^^u + di-f^d + si-f5s) (41) 



^ The relation (|38)) had also been obtained in j29|[based on a different approach. For the SU(2) version of 
Eq. (|38p . namely, (p|ui75it + dij^dlp) ~ 0, see [30|. 
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with m = {l/rriu + l/rrid + l/m^) ^. Naively, it is expected that the relation 



N 



Stt' 



GG 



N 



(42) 



permits one to recast all the unknown matrix elements in terms of the couplings g'\. However, 
it is straightforward to check that Eqs. (|33p and ()42p when combined together do not lead 
to any solution for Eq^s. As pointed out in [36|], a key to this paradox is the observation 
made in 37 1 that it is the disconnected insertion of 5£^p"'^' that is related to the insertion 



of OGG in the hadronic process 37 1. The connected insertion of 5£^p"'^' must vanish for 
on-shell amplitudes. Since a quark loop with the insertion of 5£^p"'" at zero momentum is 
the same as an insertion of OGG, the trick of rotating the 9 term into a pseudoscalar density 
does not change the nature of the calculation which must be done in the lattice or quark 
model framework [37|. Hence, if care is not taken, the use of Baluni's Lagrangian may lead 
to fake results, such as the electric dipole moments of the constituent quarks [38 1. 



4. Another piece of information on (A^|GG|A^) comes from the U(l) Goldberger-Treiman (GT) 
relation (for a review of the isosinglet GT relation, see |24i |) 



9'a 



2mi\] 



■QrioNN H 



(43) 



valid in the chiral limit, where grj^NN is a coupling constant between ?yo and the nucleon and 
the ellipsis • • • denotes possible contributions from pseudoscalar glueballs. Eq. ([I3|) is the 
analog of the familiar isovector GT relation g\ = fTr/{'2mN)gTrNN- It follows from Eq. ([33l) 
that 



N 



^GG 

Svr 



2^3 



grioNN + 



(44) 



In principle one can use this GT relation to fix the matrix element {N\GG\n) and to cross 
check the result of Eq. (fl6|) given below. In practice, this is not an easy task because of the 
mixing of t]q with i]s and and possibly the pseudoscalar glueball and the ill-determined 
couplings grj'NN and gr^NN- 

To summarize, we have listed several different methods of extracting the proton matrix elements 
of quark pseudoscalar densities and GG. It appears that only the approach based on the large- A'c 
chiral limit relation Eq. ([3]) still remains viable. Let us first quote the results from 

1 mj\[ 



and 



A^ 



^GG 
Svr 



Eu 
Ed 
Es 

N 



2 mdms{l + z + w) 

1 TTlTV 

2 mdms{l + z + w) 

1 niN 

2 mdms{l + z + w) 



{md + 2m,)g\ + mdg\ , 
-{mu + 2ms)g\ + mug\ , 
{ruu - md)gA - {rriu + md)g^^ 



(45) 



1 mN 2 Q o 
o TT"^ : ^mdm^l + z + w)gA + ms(md - mu)gA 

2 mdms[l + z + w) 13 

.li^M..^d)-2^.rn,A 



(46) 
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with z = mu/md and w = rriu/ms- Numerically, z 
pseudoscalar-nucleon coupling now has the expression 



0.5 and w = 0.025 in this work. The 



9aNN 



+ - 



h 

l + z- 



2w 



1 



z 



+ z + w 



9\ 



\ + z + w 



9a 



(47) 



In terms of Aq's and m defined in Eq. (j4ip . the matrix element {N\{as/8Tr)GG\N) and the 
pseudoscalar-nucleon coupling can be recast to more symmetric forms 



N 



8^" 



GG 



N 



771^771^(1 + Z + W 



and 

gaNN 



{V2GFf'^mN 



+ 



{V2GFf''^mN 



Au Ad As 

mu rud ms 



Au + 



(48) 




ed-(E^.)Y 



Z 



+ z + w 



Ad 



(49) 



777 



Ad + 



6 



777 
777. 



Ee. 



As 



where use has been made of Eq. (j36p . Although the relation ([3|) is subject to l/N^ and chiral 
corrections, it is pertinent to regard Eqs. (j45|) - (j49p as benchmark results. 



It is interesting to notice that the general axion-nucleon coupling has precisely the same expres- 
sion as Eq. ([^7j) or ([1^ (see 39 1 and Eq. (9) on p. 497 of [3]). However, the derivation of the axion 
coupling with the nucleons does not rely on the assumption of the large- iVc chiral limit. This has 
to do with the fact that the physical axion does not have an anomalous GG interactions 0]. 

Before proceeding further, we would like to stress that contrary to the case for scalar-nucleon 
couplings, 9crNN does receive significant contributions from light quarks even in the chiral limit (see 
Eq. (j47p or (09])). This is because the pseudoscalar boson has an admixture with vr, rfs and 770 which 
in turn couple strongly with the nucleon. The smallness of the pseudoscalar-nucleon coupling is 
compensated by the smallness of the Goldstone bosons as the ratio of niq/mTj, for example, remains 
finite in the chiral limit. 

As for the couplings 5^, we note that since there is no anomalous dimension associated with 
the axial- vector currents and A^, the non-singlet couplings g\ and g\ do not evolve with 
and hence can be determined at = from low-energy neutron and hyperon beta decays. Under 
flavor SU(3) symmetry, the non-singlet couplings are related to the SU(3) parameters F and D by 



9\ 



F + D 



9\ 



3F-D . 



(50) 



We use the updated coupling g\ = 1.2701 ± 0.0025 [jj] and the value F/D = 0.575 ± 0.016 [40|] to 
obtain 



F = 0.464 ±0.014, £»= 0.806 ± 0.008 , 5^ = 0.585 ± 0.039 . 



(51) 
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Using this value of the recent measurements of by COMPASS j41| and HERMES [4g] imply 
that 



AScoMPAss(3GeV2) 
ASH,,MEs(5GeV2) 



0.35 ±0.03 ±0.05, 

0.330 ± O.Oll(theo) ± 0.025(exp) ± 0.028(evol) . 



(52) 



Therefore, g^^ is in the vicinity of 0.34 . However, whether the flavor SU(3) symmetry can be applied 
to determine g\ still remains controversial 43l445l|. Indeed, there had been some suggestions that 
SU(3) symmetry might be badly broken and the errors could be as large as 25% |43l. l46l| . A recent 
calculation based on the cloudy bag model showed that g\ = 0.46 ± 0.05 which was reduced by 
as much as 20% below the usual phenomenological value It is clear from Eq. (j37p that a 
decrease in g\ will increase g'^ so that the magnitude of As inferred from —As = ^{g\ — g^) 
is reduced from 0.08 to 0.03 (see Table HTj). Interestingly, a very recent lattice calculation yields 
As = -0.020 ± 0.010 ± 0.004 at ^ = 2.7 GeV fl. 

There are still two remaining issues needed to be addressed. First, Eqs. (j45]) - ([49]l are derived 
based on the relation psp valid in the chiral and large- iVc limits. It is natural to ask how sizable 
the chiral and 1/Nc corrections to Eq. ([38]) are. In principle, this can be studied in the effective 
chiral Lagrangian approach that accommodates both corrections. A detailed investigation will be 
reported elsewhere. Second, since Ag's and g^ are thus far evaluated at around 4 GeV ^, 
how does one evolve them to = 0? Unfortunately, the evolution into = is beyond the 
perturbative QCD regime and hence not calculable in a perturbative way. 



Finally we come to the estimation of the coupling strength of QaNN ■ We shall assume that 



1 



for definiteness and simplicity and evaluate various parameters for two cases: (i) g\ = 0.585 and 
g^ = 0.34 and (ii) g^ = 0.46 and 5° = 0.37. ^ Note that for the neutron, {g\)n is opposite to 
{g\)p] that is, {g\)n = ~{F + D). The numerical results are summarized in Table HIl Obviously, 
the nucleon matrix elements (A^|sz75s| A^) and {N\GG\N) strongly violate isospin symmetry, as 



pointed out in 33|, [sj]. For example, it is shown in 3j] that can be an isoscalar operator only 



if the pseudoscalar operator sij^s is not an isoscalar one. It is worth mentioning that the result for 
{N\GG\N) in Table HIl will be useful if one considers the gluon field coupling to the fermionic dark 
matter through the effective operator {x'y^x)GG or to the complex scalar dark matter described 
by the interaction {x^x)GG. 



IV. SUMMARY 

The couplings between nucleons and scalar or pseudoscalar Higgs bosons are of high interest, par- 
ticularly for the estimation of cross sections between dark matter and nuclei. Such calculations rely 
on sufficient knowledge of various nucleon matrix elements of quark scalar and pseudoscalar densi- 
ties. In view of better understanding of baryon masses, quark spin components and lattice studies 
of the pion- nucleon and quark sigma terms in recent years, we update the scalar/pseudoscalar 
couplings of the nucleons in this work. 



^ The value = 0.37 is obtained from Eq. (|37)) with the input of g\ ~ 0.46 . 
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TABLE II: Summary of parameters. All are in units of MeV except for the dimensionless Ag's 
and QaNN- The values in the upper (lower) sector are obtained using = 0.586 and g'^ = 0.34 
(5,8 = 0.46 and g^ = 0.37). Note that g^ is equal to 1.2701 for the proton and -1.2701 for the 
neutron. The assumptions of three heavy quarks and = 1 have been made. The quark spin 
components Ag's and the coupling g\ are evaluated at of order 4 GeV^. 



N 


rriuEu 


nidEd 




{N\I^GG\N) 


An 


Ad 


As 


QaNN 


P 


405 


-787 


-466 


389 


0.85 


-0.42 


-0.08 


-8.2 X 10^3 


n 


-396 


796 


-75 


-2 


-0.42 


0.85 


-0.08 


1.3 X 10-3 


P 


404 


-788 


-408 


380 


0.84 


-0.44 


-0.03 


-7.8 X 10-3 


n 


-397 


795 


-17 


-11 


-0.44 


0.84 


-0.03 


1.7 X 10^3 



Early calculations usually find a larger strange quark content, characterized by the parameter y, 
and a larger value of the strange-quark sigma term as ■ They are partly due to a larger (TtiN — 55 — 60 
MeV, assumed to be the on-shell pion-nucleon scattering amplitude at the Cheng-Dashen point, 
as well as a smaller o"o — 24 MeV, extracted from the baryon masses using the leading-order 
formulae in quark masses. In this work, we take the cue from recent lattice studies that ag is much 
smaller than previously thought. For definiteness, we use as = 50 MeV as a benchmark value, 
which is consistent with the global fit of the lattice calculations. This together with cq — 36 MeV 
obtained using the baryon mass formulae computed to the second order in light quark masses leads 
to (TttTv — 40 MeV in good agreement with recent lattice determinations [19l. l23l|. Therefore, the 
strange quark content in the nucleon is y ~ 0.09 . This is roughly a factor of two smaller than the 
previous estimate of y ~ 0.2 . The relative importance of the light quarks in the effective scalar 
Higgs-nucleon coupling is shown by the ratios (T„ : cr^ : — 1 : 1.3 : 3, to be contrasted with 
1 : 1.5 : 6 in [47|] and 1 : 1.3 : 6 in [48(. (See Appendix Rl for a detailed comparison.) We find that 
the coupling g^nn — 1.1 x 10"^^ of which the strange quark contribution is 18%, a factor of two 



smaller than what is usually used in the calculations of direct detection of dark matter [47|, |48| . 
The dominant contribution (~ 70%) comes from the heavy quarks. Moreover, the ratio g^ppp/gi/mn 
is very close to unity, showing essentially no isospin violation. 

For the pseudoscalar Higgs couplings with the nucleons, several different methods for extracting 
the nucleon matrix elements of pseudoscalar quark densities and GG are reviewed. We conclude 
that the method based on the relation {Nlui-y^u + dij^d + sij5s\N) = derived in the chiral and 
large- iVc limits still remains viable. We then take two sets of axial coupling constants {g\, g^, g^) = 
(1.2701,0.585,0.34) and (1.2701,0.46,0.37) to compute the corresponding quark spin components 
Aq's. The second set of parameters is motivated by the recent observation that there is a possible 
SU(3)-breaking effect on the determination of 5^ implied by both model and lattice analyses. For 
example, a recent calculation based on the cloudy bag model finds g^ to be reduced by ~ 20% from 
the usual phenomenological value in the first set. We find that —As changes from 0.08 to 0.03 due 
to such a reduction in g^. The isospin violation in the nucleon matrix elements (A^|si75s|A^) and 
{N\GG\N) results in significantly different g^NN for N = p and n. We obtain gann/dapp ~ —0.16 
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and —0.22 for the two parameter sets, respectively. 
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Appendix A: Comparison with other works 

In the lattice and dark matter communities, it is conventional to consider the quantities defined 

by 

_ {N\mgqq\N) ""^^ , . 

rriN rriN q 

parameterizing the fractions of the nucleon mass niN- carried by the corresponding quarks and 
gluons. In terms of these parameters, the coupling g^pp (see Eq. ((28]) ) has the expression 



{V2GF)'/'mN <i fT. + ^fTa E Ch^ ■ (A2) 



9<j>pp 



A popular set of values is that used in the DarkSUSY package [4 



= 0.023 , f!f^ = 0.034 , f^] = 0.14 , f^^^ ^ 
p(") _ n mo _ n n/ii — [\ ^ a ^(") 



f}^'J = 0.019 , ffj = 0.041 , ffj = 0.14 , /^"; ^ = 0.0592 , (A3) 
and 

An(P) = Ad^") = 0.77 , Ad^^) = Au^") = -0.40 , As^^) = As^") = -0.12 . (A4) 



Another set of parameters commonly quoted in the literature is j48l | 



f}^J = 0.020 , = 0.026 , P^J = 0.118 , ^ = 0.0619 , 

fP = 0.014 , = 0.036 , = 0.118 , fP = 0.0616 , (A5) 



id 
and 



An(P) = Ad(") = 0.78 , Ad^^) = Aw^") = -0.48 , As^^) = As^") = -0.15 , (A6) 

where only the central values are listed. 
For a comparison, in this work we have 

/^P) = 0.017 , = 0.022 , 4f = 0-053 , ^ = 0.0672 , 

= 0.011 , = 0.034 , = 0.053 , ^ = 0.0668 , (A7) 

for the choice of CTg = 50 MeV, while 

An(P) = Ad^") = 0.85 , Ad^^) = Au^") = -0.42 , As^p) = As^") = -0.08 (A8) 
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for g\ = 0.585 and 



An(p) = Ad(") = 0.84 , Ad^P) = Au^") = -0.44 , As^p) = As^") 



-0.03 



(A9) 



ing s, the contributions from light quarks 



471 . |48[. In particular, the strange quark 



for g\ = 0.46 . 

It is noticed that in the scalar Higgs-nucleon coupl 
in our work are generally smaller than those used in 
contribution is nearly two to three times smaller. The heavy quarks have a more dominant effect. 
The strange-quark spin component is also found to be 2 to 5 times smaller in magnitude than those 
in HH]- 
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